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We perform a model study of deconfinement and cliiral symmetry restoration in a strong magnetic 
background. We use a Nambu-Jona Lasinio model with the Polyakov loop, taking into account a 
possible dependence of the coupling on the Polyakov loop expectation value, as suggested by the 
recent literature. Our main result is that, within this model, the deconfinement and chiral crossovers 
of QCD in strong magnetic field are entangled even at the largest value of eB considered here, namely 
eB — 30m^ (that is, B ~ 6 x 10^^ Tesla). The amount of split that we measure is, at this value of 
eB, of the order of 2%. We also study briefly the role of the 8-quark term on the entanglement of 
the two crossovers. We then compare the phase diagram of this model with previous results, as well 
as with available Lattice data. 
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I. INTRODUCTION 

The study of the Quantum Chromodynamics (QCD) 
vacuum, and of its modifications under the influence of 
external factor like temperature, baryon chemical poten- 
tial, external fields, is one of the most attractive topics of 
modern physics. One of the best strategies to overcome 
the difficulty to study chiral symmetry breaking and de- 
confinement, which share a non-perturbative origin, is 
offered by Lattice QCD simulations at zero chemical po- 
tential [iHsl- At vanishing quark chemical potential, it 
is established that two crossovers take place in a nar- 
row range of temperature; one for quark deconfinement, 
and another one for the (approximate) restoration of chi- 
ral symmetry. Besides, the use of the Schwinger-Dyson 
equations for the quark self-energy f^, 0], and the use 
of functional renormalization group [8| to the Hamilto- 
nian formulation of Yang-Mills theory in Coulomb gauge, 
are very promising QCD with one quark flavor 

at finite temperature and quark chemical potential has 
been considered in pT| within the functional renormal- 
ization group approach, combined with re-bosonization 
tecniques; in [IJI, the renormalization group fiow of the 
Polyakov-Loop potential and the flow of the chiral order 
parameter have been computed. The results of [Hj sug- 
gest that the chiral and deconfinement phase transition 
temperature agree within a few MeV for vanishing and 
(small) quark chemical potentials. Furthermore, an in- 
teresting possibility to map solutions of the Yang-Mills 
equations of motion with those of a scalar field theory 
has been proved in [l3| . 

An alternative approach to the physics of strong in- 
teractions, which is capable to capture some of the 
non-perturbative properties of the QCD vacuum, is the 
Nambu-Jona Lasinio (NJL) model [iJI, see Refs. for 
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reviews. In this model, the QCD gluon-mediated interac- 
tions are replaced by effective interactions among quarks, 
which arc built in order to respect the global symme- 
tries of QCD. Beside this, the parameters of the model 
are fixed to reproduce some phenomenological quantity 
of the QCD vacuum; therefore, it is reasonable that the 
main characteristics of its phase diagram represent, at 
least qualitatively, those of QCD. 

In recent years, the NJL model has been improved 
in order to be capable to compute quantities which 
are related to the confinement-deconfinement transi- 
tion of QCD. It is well known that color confinement 
can be described in terms of the center symmetr y o f 
the color gauge group and of the Polyakov loop [16j . 
which is an order parameter for the center symmetry 
in the pure gauge theory. In theories with dynami- 
cal fermions, the Polyakov loop is still a good indica- 
tor for the confinement-deconfinement transition, as sug- 
gested by Lattice data at zero chemical potential 

Motivated by this property, the Polyakov extended 
Nambu-Jona Lasinio model (P-NJL model) has been in- 
troduced [13 J in which the concept of statistical con- 
finement replaces that of the true confinement of QCD, 
and an effective potential describing interaction among 
the chiral condensate and the Polyakov loop is achieved 
by the coupling of quarks to a background temporal gluon 
field, and then integrating over quark fields in the parti- 
tion function. The P-NJL model, as well as its renormal- 
izable extension, namely the Polyakov extended Quark- 
Meson model (P-QM), have been studied extensively in 
many contexts jl9l - l36| . 

In a remarkable paper it has been shown by Kondo 
that it is possible to derive the effective 4-quark inter- 
action of the NJL model, starting from the QCD la- 
grangian. In his derivation, Kondo has shown explicitly 
that the NJL vertex has a non-local structure, that is, it 
is momentum-dependent; besides, the vertex acquires a 
non-trivial dependence on the phase of the Polyakov loop. 
This idea has been implemented within the P-NJL model 
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in [3l|; the modified model has then been named EPNJL, 
and the Polyakov-loop-dependent vertex has been called 
entanglement vertex. We will make use of this nomencla- 
ture in the present article. Before going ahead, it worths 
to notice that a low-energy limit of QCD, leading to the 
non-local Nambu- Jona Lasinio model studied in [38| , has 
been discussed independently in [39| . In this reference, 
the low-energy limit of the gluon propagator leads to a 
relation among the NJL coupling constant and the string 
tension. 

In this article, we report on our study of deconfine- 
mcnt and chiral symmetry restoration at finite temper- 
ature in a strong magnetic background. This study is 
motivated by several reasons. Firstly, it is extremely in- 
teresting to understand how an external field can mod- 
ify the main characteristics of confinement and sponta- 
neous chiral symmetry breaking. Lattice studies on QCD 
in magnetic (as well as chromo-magnetic) backgrounds 
can be found in [40l - !42| . Studies of QCD in magnetic 
fields, and of QCD-like theories as well, can be found 
in Refs. [43l - l50l |. Besides, strong magnetic fields might 
be produced in non-central heavy ion collisions jsil [52 | . 
More concretely, at the center-of-mass energy reachable 
at LHC, ^/snn ~ 4.5 TeV, the magnetic field can be as 
large asQ ~ 15m^ according to It has been ar- 
gued that in these conditions, the sphalcron transitions 
of finite temperature QCD, give rise to Chiral Magnetic 
Effect (CME) mil. 

The novelty of this study is the use of the EPNJL 
model in our calculations, in the one-loop approximation. 
In comparison with the original PNJL model of [l^ , the 
EPNJL model has two additional parameters. However, 
they are fixed from the QCD thermodynamics at zero 
magnetic field, as we will discuss in more detail later. 
Therefore, the results at a finite value of the magnetic 
field strength have to be considered as predictions of 
the model. Our main result is that the entanglement 
of the NJL coupling constant and the Polyakov loop, 
might affect crucially the phase diagram in the tempera- 
ture/magnetic field strength plane. Previous model stud- 
ies |43l - l45| have revealed that both the deconfinemcnt 
temperature, T^, and the chiral symmetry restoration 
temperature, T^, are enhanced by a magnetic field, in 
agreement with the Lattice data of Ref. [40| ■ The model 
results are in slight disagreement with the Lattice, in the 
sense that the former predict a considerable split of Tl 
and as the strength of the magnetic field is increased. 
Within the EPNJL model, we can anticipate one of the 
results, that is, the split among Tl and might be con- 
siderably reduced even at large values of the magnetic 
field strength. In particular, using the values of the pa- 
rameters of [3l|, which arise from a best- fit of Lattice 
data at zero and imaginary chemical potential and which 



In this article, wc measure eB in units of the vacuum squared 
pion mass m^; then, eB = rn^ corresponds to _B 2 X 10"'^'* 
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are appropriate for our study, we find a split of the or- 
der of 2% at the largest value of eB considered, namely 
eB = 30toJ. 

In (45j a similar computation within a model without 
entanglement, but with an 8-quark interaction added, has 
been performed. Following the nomenclature of [3l|, we 
call the latter model, PNJLg model. We find the com- 
parison among the EPNJL and the PNJLg models very 
instructive. As a matter of fact, the parameters in the 
two models are chosen to reproduce the QCD thermody- 
namics at zero and imaginary chemical potential [28l . [31| . 
Therefore, both of them are capable to describe QCD in 
the same regime. It is interesting that, because of the 
different interactions content, the two models predict a 
slight different behavior of hot quark matter in strong 
magnetic field. In the next future, the comparison with 
refined Lattice data can enlighten on which of the two 
models is a more faithful description of QCD. 

The paper is organized as follows: in Section II we sum- 
marize the formalism: we derive the quark propagator 
and the equation for the chiral condensates in magnetic 
field, within the EPNJL model, in the one-loop approx- 
imation; then, we compute the thermodynamic poten- 
tial, in Section III, we collect our results for the chiral 
condensate and the expectation value of the Polyakov 
loop, for several values of the magnetic field strength. In 
Section IV, we briefiy investigate on the effect of the 8- 
quark interaction in the EPNJL model in magnetic field. 
In Section V, we draw the phase diagram of the EPNJL 
model in magnetic field, and make a comparison with our 
previous result [45| . Finally, in Section VI we draw our 
conclusions, and briefly comment on possible extensions 
and prosecutions of our study. 

We use natural units throughout this paper, ft = c = 
ks = 1, and work in Euclidean space-time i?"* = (5V 
, where V is the volume and (3 = \/T with T corre- 
sponding to the temperature of the system. Moreover, 
we take a non-zero current quark mass. In this case, 
both deconfinemcnt and chiral symmetry breaking are 
crossovers; however, we sometimes will make use of the 
term "phase transition" to describe them, for stylistic 
reasons. It should be clear from the context that our 
"phase transitions" are meant to be crossovers unless 
stated differently. 



II. THE MODEL WITH ENTANGLEMENT 
VERTEX 

We consider a model in which quark interaction is de- 
scribed by the following lagrangian density: 

C^^ {i-1^D^ - m) + £/ ; (1) 

here ^ is the quark Dirac spinor in the fundamental repre- 
sentation of the fiavor SU{2) and the color group; r cor- 
respond to the Pauli matrices in fiavor space. A sum over 
color and flavor is understood. The covariant derivative 
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embeds the QED coupling of the quarks with the exter- 
nal magnetic field, as well as the QCD coupling with the 
background gluon field which is related to the Polyakov 
loop, see below. Furthermore, we have defined 



Cj=G 



(2) 



This interaction term is invariant under SU{2)y <Si 
SU{2)a (E) U{1)v- In the chiral limit, this is the sym- 
metry group of the action as well, if no magnetic field 
is applied. However, this group is broken explicitly to 
U{l)y (g) U{1)\ (g) U{l)v if the magnetic field is coupled 
to the quarks, because of the different electric charge of 
u and d quarks. Here, the superscript 3 in the V and 
A groups denotes the transformations generated by T3, 
''375 respectively. Therefore, the chiral group in presence 
of a magnetic field is U (g) U{1)\. This group is then 
explicitly broken by the quark mass term to 

We are interested to the interplay among chiral sym- 
metry restoration and deconfinement in a strong mag- 
netic field. To compute a temperature for the decon- 
finement crossover, we use the expectation value of the 
Polyakov loop, that we denote by L. In order to com- 
pute L we introduce a static, homogeneous and Euclidean 
background temporal gluon field, Aq = iA^ = iXaA^, 
coupled minimally to the quarks via the QCD covariant 
derivative O. Then 



study, we follow a more phenomenological approach to 
the problem, using the ansatz introduced in j31| . that is 



G = g [1 - aiLL - a2{L^ + L^)] 



(5) 



and we take L = L from the beginning. The functional 
form in the above equation is constrained by C and ex- 
tended Z3 symmetry. We refer to [U for a more detailed 
discussion. The numerical values of ai and 0,2 have been 
fixed in [3l| by a best fit of the available Lattice data at 
zero and imaginary chemical potential of Ref. [54l | , which 
have been confirmed recently in [ssf . In particular, the 
fitted data are the critical temperature at zero chemi- 
cal potential, and the dependence of the Roberge- Weiss 
endpoint on the bare quark mass. 

The values ai = a2 = a = 0.2 ± 0.05 have been ob- 
tained in [31[ using a hard cutoff regularization scheme. 
We will focus mainly on the case a = 0.2 as in 3l|. We 
have verified that in our regularization scheme, our re- 
sults are in quantitative agreement with those of jsH ]. 
when we take the parameter Tq = 190 MeV in the 
Polyakov loop effective potential in agreement with [3l[ , 
see below. Then, we will study how the results change 
when we vary a. 



A. Quark propagator and chiral condensate in 
magnetic field 



L= -Tr,exp(2;3Aa^4) 



(3) 



where j3 = \/T. In the Polyakov gauge, which is con- 
venient for this study, ^0 = iX?,4> + iAs^^; moreover, 
we work at zero quark chemical potential, therefore we 
can take L = from the beginning, which implies 
^4 = 0. This choice is also motivated by the study 
of [131, where it is shown that the paramagnetic con- 
tribution of the quarks to the thermodynamic potential 
induces the breaking of the Z3 symmetry, favoring the 
configurations with a real-valued Polyakov loop (see Sec- 
tion III.C of |44| for an excellent discussion of this point). 
We have then[ll, 



L = 



1 + 2 cos(/30) 



(4) 



As already discussed in the Introduction, it has been 
shown that it is possible to derive the effective 4-quark 
interaction ^ starting from the QCD lagrangian j37l |. 
In [131 it has been shown that the NJL vertex has a non- 
local structure, that is, it is momentum-dependent. An 
analogous conclusion is achieved in [s^ . More important 
for our study, the NJL vertex acquires a non-trivial de- 
pendence on the phase of the Polyakov loop. Therefore, 
in the model we consider here, it is important to keep 
into account this dependence. The exact dependence of 
G on L has not yet been computed; it is possible that 
it will be determined in the next future by means of the 
functional renormalization group approach [37j . In our 



We work in the mean field approximation throughout 
this paper, neglecting pseudoscalar condensates; more- 
over, we make the assumption that condensation takes 
place only in the fiavor channels tq and T3. The mean 
field interaction term Eq. Q can be cast in the form 



C = -2GE {uu + dd) - GY? 



(6) 



where E = — {uu + dd) . 

In a magnetic field, the chiral condensates of u and d 
quarks have to be different, because the electric charges 
of these quarks are different. Even if the one- loop quark 
self-energies in Eq. ([6]) depend on the sum of the two 
condensates, being therefore flavor independent, it is 
straightforward to show that the two condensates turn 
out to be different, by taking the trace of the propagator 
of the two quarks. The interaction ^ is diagonal in fla- 
vor space, therefore we can focus on the propagator of a 
single flavor /. 

To write the quark propagator we use the Ritus 
method [56j . which allows to expand the propagator on 
the complete and orthonormal set made of the eigenfunc- 
tions of a charged fermion in a homogeneous and static 
magnetic field. This is a well known procedure, discussed 
many times in the literature, see for example (57l - [60| : 
therefore it is enough to quote the final result. 



k=0 



dpodp2dp3 
(27r)4 



1 



Ep{v) , 
(7) 
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where Ep{x) corresponds to the eigenfunction of a 
charged fermion in magnetic field, and Ep{x) = 
7o(i?p(x))^7o. In the above equation, 



P = {p^+iAi,Q,Qj2k\QfeB\,p,) , 



(8) 



where k = 0,1,2,... labels the fc*^ Landau level, and 
Q = sign{Qf), with Qf denoting the charge of the fla- 
vor /; Afe is a projector in Dirac space which keeps into 
account the degeneracy of the Landau levels; it is given 
by 

Ak = 4o [V+Sq^+i + V-Sq,-i] + (1 - dko)I , (9) 

where V± are spin projectors and / is the identity matrix 
in Dirac spinor indices. The propagator in Eq. ([7]) has 
a non-trivial color structure, due to the coupling to the 
background gauge field, see Eq. ([5]). 

The trace of the /-quark propagator is minus the chiral 
condensate (//), with / = u,d. Taking the trace in 
coordinate and internal space, it is easy to show that the 
following equation holds: 



iff) = -N. 



\QfeB\ 
2tt 



■f2pJ^^CiL,L,T\p.,k) 



(10) 



Here, 



C{L,L,T\p„k) = UA~2U{L,L,T\p,,k) , (11) 



and M denotes the statistically confining Fermi distribu- 
tion function. 



C{L,L,T\p,,k) 



1 + 2Le^"^ + Le^l^'^f 



where 



with 



j} =pl + 2\QfeB\k + M] 



Mu = Md = mo + 2GS . 



(12) 
(13) 
(14) 



The first and the second addenda in the r.h.s. of Eq. 
correspond to the vacuum and the thermal fluctuations 
contribution to the chiral condensate, respectively. The 
coefficient (3^ — 2 — Sko keeps into account the degen- 
eracy of the Landau levels. The vacuum contribution is 
ultraviolet divergent. In order to regularize it, we adopt 
a smooth regulator U\, which is more suitable, from the 
numerical point of view, in our model calculation with 
respect to the hard-cutoff which is used in analogous cal- 
culations without magnetic field. In this article we chose 



Ua 



A 



2N 



A2^ + (p2+2|Q/eB|fc) 



N 



(15) 



To be specific, we consider here the case N = 5, for 
numerical convenience. A larger value of N increases the 
cutoff artifacts, as already discussed in detail in [H, IH, 
|47| , but leaves the qualitative picture unchanged; on the 
other hand, a smaller value of N, namely iV < 3, makes 
not possible the fit of the pion decay constant and of the 
chiral condensate in the vacuum. The (more usual) 3- 
momcntum cutoff rcgularization scheme is recovered in 
the limit N — > oo; we notice that, even if the choice 
N = 5 may seem arbitrary to some extent, it is not more 
arbitrary than the choice of the hard cutoff, that is, of a 
regularization scheme. 

By virtue of Eq. ([TU| it is easy to argue that the con- 
densates of u and d quarks are different in magnetic field. 
As a matter of fact, the different value of the charges 
makes the r.h.s. of the above equation flavor-dependent, 
even if S is flavor-singlet. Once we compute the mean 
field value of E by minimization of the thermodynamic 
potential (see the next Section), we will use Eq. ([TU)) to 
evaluate the chiral condensates for u and d quarks in 
magnetic field. 



B. Thermodynamic potential 

In the one-loop approximation, the thermodynamic 
potential is given by 

n=U{L,L,T) + GJ:^ - ^Trlog {f3S-') , (16) 

where the trace is over color, flavor, Dirac and space- 
time indices and the propagator for each flavor is given 
by Eq. ([7]). It is straightforward to derive the final result, 
namely 



n=U{L,L,T) + UM 

^ f+^dp 
' 1^ -^T^l^^k / -:T-G{L,L,T\p,,k) . 

(17) 

In the above equation, the first line corresponds to the 
classical contribution; we have defined 



(18) 



the second line corresponds to the sum of vacuum and 
thermal one loop contributions, respectively, and arise af- 
ter integration over the fermion fluctuations in the func- 
tional integral. We have defined 



g(L, L, T\p,, k) = N,UA{p)iOf + - log J- , (19) 



with 



J-(L, Z, T\p,, fc) = 1 + 3Le-'^"^ + 3Le-2'5"^' + e'^^^^ . 

(20) 
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The potential term U[L,L,T] in Eq. ^ is built 
by hand in order to reproduce the pure gluonic lat- 
tice data [i^ [131 ■ We adopt the following logarithmic 
form [ii, 



U[L,L,T] = - 



a{T) 



LL 



+ b{T) ln[l - 6LL + 4(1^ + L^) - 2,{LLf] | , 

(21) 

with three model parameters (one of four is constrained 
by the Stefan-Boltzmann limit), 



a{T) = ao + ai ( ^ j + aa 



(22) 



The standard choice of the parameters reads ag = 3.51, 
ai = —2.47; a2 ~ 15.2 and 63 = —1.75. The parameter 
To in Eq. (j2ip sets the deconfinement scale in the pure 
gauge theory. In absence of dynamical fermions one has 
To = 270 MeV. However, dynamical fermions induce a 
dependence of this parameter on the number of active 
flavors [3^ . For the case of two light flavors to which we 
are interested here, we take Tq = 212 MeV as in [45|. 

As a final remark, it is easy to show that summing 
Eq. PU)) for u and d quarks, one reproduces the equa- 
tion satisfied by E, the latter being obtained from the 
stationarity condition dfl/dH = 0. 
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FIG. 1. Upper panel: Chiral condensate E, measured in units 
of the chiral condensate at zero temperature and zero field, 
Eo, as a function of temperature, for several values of the 
magnetic field strength. Lower panel: Expectation value of 
the Polyakov loop, L, as a function of temperature, for several 
values of the magnetic field strength. The data are obtained 
for a = 0.2. In the figures, the magnetic fields are measured 
in units of m?.. 



III. DECONFINEMENT AND CHIRAL 
SYMMETRY RESTORATION 

In this Section wc summarize our numerical results. In 
this study we fix the values of A and mo to reproduce 
the values of /^r and TOtt in the vacuum, as well as the 
numerical value of the light quarks chiral condensate. Wc 
have then A = 626.76 MeV, g = 2.02/A^ The pion 
mass in the vacuum, TOtt = 139 MeV, is used to fix the 
numerical value of the bare quark mass via the GMOR 
relation /^m^ = —2mo{uu); this gives mg = 5 MeV 
when we take (uu) = (—253 MeV)^. Finally, we take 
To = 212 MeV in the Polyakov loop effective potential as 
in [4^, unless stated differently The values of E and L 
are computed by the minimization of the thermodynamic 
potential. 

A. The case a-i = 02 = 0.2 

In Fig. [T] we plot our data for the chiral condensate 
E, measured in units of the condensate at T = eB = 0, 
that is Eq = 2 X (—253 MeV)^, and the expectation value 
of the Polyakov loop as a function of temperature, com- 
puted for several values of the magnetic field strength. 



They are computed from the minimization of the ther- 
modynamic potential in Eq. p7p . In the numerical com- 
putation, we have chosen ai = 02 = a = 0.2 as in [3l| . 

The results shown in Fig. [1] arc interesting for sev- 
eral reasons. Firstly, if we identify the deconfinement 
crossover with the temperature T^ at which dL/dT is 
maximum, and the chiral crossover with the temperature 
at which |dE/dT| is maximum, we observe that the 
two temperatures are very close also in a strong magnetic 
field. For concreteness, at eB = we find Tp = = 
185.5 MeV. Besides, at eB = 19ml we find Tp = 199 
MeV and T^ = 201 MeV. Hence at the strongest value 
of eB considered here, namely eB = 19m^, the entan- 
glement vertex makes the split of the two crossovers 
of 1.5%. From the model point of view, it is easy 
to understand why deconfinement and chiral symmetry 
restoration are entangled also in strong magnetic field. 
As a matter of fact, using the data shown in Fig. [1] it 
is possible to show that the NJL coupling constant in 
the pseudo-critical region in this model decreases of the 
15% as a consequence of the deconfinement crossover. 
Therefore, the strength of the interaction responsible for 
the spontaneous chiral symmetry breaking is strongly 
affected by the deconfinement. with the obvious conse- 
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quence that the numerical value of the chiral condensate 
drops down and the chiral crossover takes place. We 
have verified that the picture remains qualitatively and 
quantitatively unchanged if we perform a calculation at 
eB = 30ml. In this case, we find Tl = 224 MeV and 
= 225 MeV. 

Before going ahead, it worths to notice that we have 
also considered the case Tq = 190 MeV, which corre- 
sponds to the value considered in [3l|. In this case, for 
eB = we find = Tl = 175 MeV, in excellent agree- 
ment with [3l|. This results is comforting, because it 
shows that even using a different regularization scheme, 
the UV-regulator does not affect the physical predic- 
tions at eB = 0. Moreover, at eB/ml = 30, we find 
T^^Tl = 215 MeV. 

This result can be compared with our previous cal- 
culation [4^ in which we did not include the Polyakov 
loop dependence of the NJL coupling constant. In [i^ 
we worked in the chiral limit and we observed that the 
Polyakov loop crossover in the PNJL model is almost in- 
sensitive to the magnetic field; on the other hand, the 
chiral phase transition temperature was found to be very 
sensitive to the strength of the applied magnetic field, in 
agreement with the well known magnetic catalysis sce- 
nario [4^. This model prediction has been confirmed 
within the Polyakov extended quark- meson model in |44| , 
when the contribution from the vacuum fcrmion fiuctua- 
tions to the energy density is kcpt into account 0; we then 
obtained a similar result in [45| , in which we turned from 
the chiral to the physical limit at which uItt ~ 139 MeV, 
and introduced the 8-quark term as well (PNJLg model, 
according to the nomenclature of [3lj|). The comparison 
with the results of the PNJLg model of [i^ is interest- 
ing because the model considered there, was tuned in 
order to reproduce the Lattice data at zero and imag- 
inary chemical potential [l^, like the model wc use in 
this study. Therefore, they share the property to de- 
scribe the QCD thermodynamics at zero and imaginary 
chemical potential; it is therefore instructive to compare 
their predictions at finite eB. 

For concreteness, in [4^ we found Tp = 185 MeV and 
— 208 MeV at eB = 19mJ, corresponding to a split 
of w 12%. On the other hand, in the present calculation 
we measure a split of ~ 1.5% at the largest value of eB 
considered. Therefore, the results of the two models are 
in slight quantitative disagreement; this disagreement is 
then reflected in a slightly different phase diagram. We 
will draw the phase diagram of the two models in a next 
Section; however, since now it is easy to understand what 
the main difference consists in: the PNJLg model pre- 
dicts some window in the eB — T plane in which chiral 



^ If the vacuum corrections are neglected, the deconfinement and 
chiral crossovers are found to be coincident even in very strong 
magnetic fields |44|| . but the critical temperature decreases as a 
function of eB; this scenario is very interesting theoretically, but 
it seems it is excluded from the recent Lattice simulations |40| . 
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FIG. 2. Upper panel: Chiral condensates of u and d quarks 
as functions of temperatures in the pseudo-critical region, at 
eB = 15ml and eB = 30m^. Condensates are measured in 
units of their value at zero magnetic field and zero tempera- 
ture, namely ao = (—253 MeV)^. Lower panel: Polyakov loop 
expectation value as a function of temperature, at eB = 15m^ 
and eB = 30ml. Data correspond to q = 0.2. 



symmetry is still broken by a chiral condensate, but de- 
confinement already took place. In the case of the EP- 
NJL model, this window is shrunk to a very small one, 
because of the entanglement of the two crossovers at fi- 
nite eB. On the other hand, it is worth to stress that 
the two models share an important qualitative feature: 
both chiral restoration and deconfinement temperatures 
are enhanced by a strong magnetic field, in qualitative 
agreement with the existing Lattice data [40| . 

A further interesting feature of our data is that in the 
low temperature region, the Polyakov loop is suppressed 
as we increase the strength of the magnetic field; on the 
other hand, in the deconfinement phase, L is enhanced 
by the magnetic field. This result was also found in our 
previous studies [i^, with a fixed coupling constant; 
more remarkably, this result is in agreement with the 
Lattice data [i^l • At the moment we do not have physical 
arguments to explain this result, but we agree with the 
authors of Ref. [40| that this aspect should be considered 
in more detail from the theoretical point of view. 

For completeness, we plot in Fig [5] the chiral conden- 
sates of u and d quarks as a function of temperature, 
at eB = 15m J and eB = 30to^. The condensates are 
measured in units of their value at zero magnetic field 
and zero temperature, namely ctq = (uu) = (dd) = 
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(—253 McV)'^. They are computed by a two-step pro- 
cedure: firstly we find the values of S and L that min- 
imize the thermodynamic potential; then, we make use 
of Eq. (jlOp to compute the expectation values of uu and 
dd in magnetic field. If we measure the strength of the 
crossover by the value of the peak of |dE/c?T|, it is ob- 
vious from the Figure that the chiral crossover becomes 
stronger and stronger as the strength of the magnetic 
field is increased, in agreement with [40l |. 



scheme is used. Therefore, we can expect that a slightly 
different value of a is needed in our case. However, it 
is useful to observe that in the limit of zero field, our 
regulator is very similar, quantitatively speaking, to the 
hard-cutoff of [31|, because of the strong power-law de- 
cay of Ua as IpI « A. Therefore, the possible numerical 
deviation in the parameters is expected to be very tiny, 
if any. We have a further corroboration we are on target 
with parameters: if we chose Tq = 190 MeV as in [3l| . 
we find = = 175 MeV, in quantitative agreement 
with [3l| : this convinces that our ultraviolet re gula tor 
does not lead to quantitative discrepancy with j3l| at 
zero magnetic field. 

We have investigated the dependence of the pseudo- 
critical temperatures as a function of a, at eB = and 
eB = 19™^. The results of our computations arc col- 
lected in the upper panel of Fig. [31 At a = 0, which cor- 
responds to the original PNJL model, we measure a split 
of the two critical temperatures of w 18% even at eB = 0. 
This is well known result in the PNJL literature [H, [13] , 
and to overcome this problem, several solutions have been 
suggested, like the use of the 8-quark interaction [23, [3l| . 
However, the discrepancy of the critical temperatures is 
enhanced as the value of eB is increased. At eB = 19m^ 
we find a split of ~ 30%. As we increase a, the two tem- 
peratures get closer rapidly, both at zero and at non-zero 
value of the magnetic field strength. At a = 0.2, which 
is the result that we have discussed previously and is the 
value quoted in [3lj |. the two temperatures are almost 
coincident. 

An interesting point that we have found is that in- 
creasing further the value of a, the crossovers become 
stronger and stronger; there exists a critical value of a at 
which the crossover is replaced by a sudden jump in the 
chiral condensate and in the Polyakov loop expectation 
value. This value of a is denoted by a green square in 
Fig. [21 and it is ei?— dependent. For concreteness, in the 
lower panel of Fig. [3] we plot our data for the half of the 
chiral condensate and the Polyakov loop at a = 0.275, 
which should be compared with the data in Fig. [1] where 
a = 0.2. 



IV. THE EFFECT OF THE 8-QUARK 
INTERACTION 



B. Varying a 

Since our study is based on a simple model, and we 
can not compute the exact values of the parameter a 
in Eq. ([S]) starting from first principles, we find very in- 
structive to measure the robustness of our results as we 
vary the numerical value of this parameter. We empha- 
size that in [3T| an estimate for the numerical value of 
a is achieved via the best fit of the model data with 
the Lattice data at zero and imaginary chemical poten- 
tial. However, in that reference a different rcgularization 



We have briefly investigated on the role of other inter- 
actions on the entanglement of deconfincment and chiral 
symmetry restoration crossovers. We report here the re- 
sults related to the effect of the 8-quark term, which has 
been extensively studied in [6l| in relation to the vacuum 
stability of the NJL model, to the hadron properties and 
to the critical temperature at zero baryon density; be- 
sides, it has been studied in the context of the P-NJL 
model in [2^, [U, [4^ . To this end, we add to the interac- 
tion lagrangian in Eq. ([2]), the following term 



(23) 
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In principle, we expect that the constant Gs acquires a 
dependence on the Polyakov loop as well; however, for 
simplicity we neglect this dependence here, since we are 
only interested to understand the qualitative effect of the 
interaction ([2^. We remark that the model studied in 
this Section is different from the P-NJLg model defined 
in [3l|; in the latter, the dependence of G on L is ne- 
glected. 

At the one-loop order we have 

C = -2uu [GE + 2G8S^] - 2dd [GE + 2GsX'^] 

-SGgE'^ - GE^ , (24) 

instead of Eq. ^ . The thermodynamic potential is still 
given by Eq. (|f 7p . with 



and 



Um = SGgE^ + GsE^ , (25) 



M„ = Md = mo-f2(GE + 2G8E3) , (26) 



instead of Eq. . The parameters are fixed as in [43] , 
that is A = 589 MeV, g = 5 x 10"^ MeV-^ Gg = 6 x 
10-22 MeV"** and mo = 5.6 MeV. 

In Fig. 0] we plot our data for the chiral condensate 
and the expectation value of the Polyakov loop, as a 
function of temperature, for several values of the mag- 
netic field strength and for a = 0.1. We find that even 
at this smaller value of a, the two crossovers are entan- 
gled at large magnetic fields. On the other hand, we find 
that taking a running NJL coupling makes the crossovers 
sharper, and eventually we measure a discontinuity in the 
order parameters for a large magnetic field strength. 

In Fig.[S]we plot the critical temperatures as a function 
of a, for eB = and eB = lOm^. At eB ^ the two 
crossovers are entangled already at a = 0; therefore, the 
effect of the entanglement vertex is to make the crossover 
sharper. Eventually, at the critical value a = 0.125 the 
crossover is replaced by a sudden jump of the expecta- 
tion values, analogously to the results that we find in 
the model without 8-quark term. As we increase eB, the 
deconfinement and chiral crossovers are spli t at a = 0, 
in agreement with our previous findings [45| . However, 
the split in this case is more modest than the split that 
we measure in the model with Gg = 0. In that case, at 
eB = IQm^ we find a split of « 30%, to be compared 
with the one that we read from Fig. [5] for the model with 
8-quark interaction, namely ~ 12%. Therefore, our con- 
clusion is that the 8-quark interaction helps to keep the 
two crossovers close in a strong magnetic field, even if 
it is not enough and the entanglement vertex has to be 
included, as it is clear from Fig. [S] 



V. THE PHASE DIAGRAM 

In Fig.|n]we collect our data on the pseudo-critical tem- 
peratures for deconfinement and chiral symmetry restora- 
tion, in the form of a phase diagram in the eB — T plane. 



1.5 



1.0 



0.5 



0.0^ 
160 



I 1 
i ! 



170 



180 190 
T (MeV) 



200 



210 




160 170 180 190 200 210 
T (MeV) 



FIG. 4. Upper panel: Chiral condensate E, measured in units 
of the chiral condensate at zero temperature and zero field, 
Eo, as a function of temperature, for several values of the 
magnetic field strength. Lower panel: Expectation value of 
the Polyakov loop P as a function of temperature, for several 
values of the magnetic field strength. The results are obtained 
within the model with 8-quark interaction and a = 0.1. In 
the figures, the magnetic fields are measured in units of m^. 
Lines with the same dashing correspond to the same value of 
the magnetic field strength. 



240 r 
220 



with 8-quark interaction 



> 



Tx,eB=19 



200 



180 



160L 



TP, e8= 19 



eB=0. 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 



FIG. 5. Pseudo-critical temperatures as functions of the pa- 
rameter a in the NJL coupling constant. Dashed lines cor- 
respond to the Polyakov loop crossover, for eB = and 
eB — 19. Dot-dashed lines correspond to the chiral crossover. 
The grey squares denote the value of a at which the crossover 
becomes a first-order phase transition. 



9 



In the upper panel we show the results obtained within 
the EPNJL model; in the lower panel, we plot the results 
of the PNJLg model, that are obtained using the fitting 
functions computed in [45| . In the figure, the magnetic 
field is measured in units of m^; temperature is mea- 
sured in units of the deconfinement pseudo-critical tem- 
perature at zero magnetic field, namely Tb=o = 185.5 
MeV for the EPNJL model, and Tb=o = 175 MeV for 
the PNJLg model. For any value of eB, we identify the 
pseudo-critical temperature with the peak of the effective 
susceptibility. 

It should be kept in mind, however, that the definition 
of a pseudo-critical temperature in this case is not unique, 
because of the crossover nature of the phenomena that 
we describe. Other satisfactory definitions include the 
temperature at which the order parameter reaches one 
half of its asymptotic value (which corresponds to the 
T limit for the chiral condensate, and to the T — )• 
+00 for the Polyakov loop) , and the position of the peak 
in the true susceptibilities. The expectation is that the 
critical temperatures computed in these different ways 
differ from each other only of few percent. This can be 
confirmed concretely using the data in Fig. [5] at eB — 
30toJ. Using the peak of the effective susceptibility we 
find = 225 MeV and Tl = 224 MeV; on the other 
hand, using the half-value criterion, we find = 227 
MeV and Tl = 222 MeV, in very good agreement with 
the previous estimate. Therefore, the qualitative picture 
that we derive within our simple calculational scheme, 
namely the entanglement of the two crossovers in a strong 
magnetic field, should not be affected by using different 
definitions of the critical temperatures. 

Firstly we focus on the phase diagram of the EPNJL 
model, which is one of the novelties of our study. In 
the upper panel of Fig. [SJ the dashed and dot-dashed 
lines correspond to the deconfinement and chiral sym- 
metry restoration pseudo-critical temperatures, respec- 
tively. We have fit our data using the following functional 
form: 



Tr 



= l + A 



\eB\ 
T2 



(27) 



where the subscripts Xi ^ correspond to the (approx- 
imate) chiral restoration and deconfinement tempera- 
tures, respectively. The functional form is taken in anal- 
ogy to the one used in [i^l . Numerical values of the best- 
fit coefficients are collected in Table ID As a consequence 
of the entanglement, the two crossovers stay closed also in 
very strong magnetic field, as we have already discussed 
in the previous Sections. The grey region in the Figure 
denotes a phase in which quark matter is (statistically) 
deconfined, but chiral symmetry is still broken. Accord- 
ing to [gI, H^, we can call this phase as Constituent 
Quark Phase (CQP). 

On the lower panel of Fig. [5] we have drawn the 
phase diagram for the PNJLg model; it is obtained us- 
ing Eq. ((27| with the coefficients computed in Ref. j45l |. 
and collected in Table HI The critical temperature at 
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FIG. 6. Upper panel: Phase diagram in the eB — T plane 
for the EPNJL model. Temperatures on the vertical axis are 
measured in units of the pseudo-critical temperature for de- 
confinement at eB — 0, namely Tc = 185.5 MeV. Lower panel: 
Phase diagram in the eB — T plane for the PNJLg model. 
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the pseudo-critical temperature for deconfinement at eB = 0, 
namely Tc = 175 MeV. In both the phase diagrams, T^, Tl 
correspond to the chiral and deconfinement pseudo-critical 
temperatures, respectively. The grey shaded region denotes 
the portion of phase diagram in which hot quark matter is 
deconfined and chiral symmetry is still broken spontaneously. 



TABLE I. Coefficients of the fit function defined in Eq. (gZ]) 
for the two models. 
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(MeV) 
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zero field in this case is a little smaller than the same 
temperature for the EPNJL model, because the param- 
eters needed to fit the vacuum properties are different. 
However, a comparison among the two phase diagrams 
is still instructive. The most astonishing feature of the 
phase diagram of the PNJLg model is the entity of the 
split among the deconfinement and the chiral restoration 
crossover. The difference with the result of the EPNJL 
model is that in the former, the entanglement with the 
Polyakov loop is neglected in the NJL coupling constant. 



10 



As we have already mentioned in the previous Section, 
the maximum amount of split that we find within the 
EPNJL model, at the largest value of magnetic field con- 
sidered here, is of the order of 2%; this number has to 
be compared with the split at eB = 20?7i^ in the PNJLg 
model, namely « 12%. The larger split causes a consid- 
erable portion of the phase diagram to be occupied by 
the CQP. 

Keeping into account the crudeness of the two mod- 
els, it is fair to say that they share one important fea- 
ture: both deconfinement and chiral symmetry restora- 
tion temperature are increased by a strong magnetic 
field. They disagree quantitatively, in the sense that the 
amount of split in the EPNJL model is more modest 
in comparison with that obtained in the PNJLg model. 
Before closing, we notice that a similar phase diagram 
calculation has been performed within the Quark-Meson 
(QM in the following) model in [i^]- In this reference, the 
magnetic field has the same effect that we measure in the 
PNJLg model, that is, to split of deconfinement and chi- 
ral symmetry restoration crossovers. Besides, both 
and Ti are enhanced by a magnetic field, in excellent 
agreement with the behavior that we find within the sev- 
eral NJL-like models that we have used in this study. 

On the Lattice side, the most recent data about this 
kind of study are those of Rcf. [13] . In these data, the 
split among the two crossovers is absent. Then the Lat- 
tice data seem to point towards the phase diagram of the 
EPNJL model. On the other hand, the results of [1^ 
might not be definitive, in the sense that both the lat- 
tice size might be enlarged, and the pion mass could be 
lowered to its physical value in the vacuum. Therefore, 
it will be interesting to compare our results with more 
refined data. 



VI. CONCLUSIONS 

In this article, we have studied chiral symmetry 
restoration and deconfinement in a strong magnetic back- 
ground, using an effective model of QCD. In particular, 
we have reported our results about the effect of the entan- 
gled vertex on the phase diagram. Our main result is that 
the entanglement reduces considerably the split among 
the deconfinement and the chiral symmetry crossovers 
studied in j43l - |45| . as expected. 

We have also studied the effect of the 8-quark term on 
the split. Our results suggest that the 8-quark interaction 
helps the two crossovers to be close. Furthermore, we 
have shown that the crossovers become sharper and then 
they are replaced by a sudden jump of the expectation 
values, as the value of a in the entanglement vertex is 
larger than a critical value. 

We have then compared our results with those of other 
model calculations, namely the PNJL model [ij . the 
Quark-Meson model ^ and the PNJLg model 0. The 
most striking similarity among the several models is that 
they all support the scenario in which chiral symme- 



try restoration and deconfinement temperatures are en- 
hanced by a strong magnetic field. The models differ 
quantitatively for the amount of split measured (very few 
percent for the EPNJL model, and of the order of 10% 
for the other models for eB « 20to^). 

Furthermore, we have compared our results with those 
obtained on the Lattice . In [4^ , the largest value of 
magnetic field considered is eB « 0.75 GeV^, which cor- 
responds to eB/rn^ « 38. The Lattice data seem to point 
towards the phase diagram of the EPNJL model. On the 
other hand, the results of [40| might not be definitive: the 
lattice size might be enlarged, the lattice spacing could 
be taken smaller (in [i^l the lattice spacing is a = 0.3 
fm), and the pion mass could be lowered to its physical 
value in the vacuum. As a consequence, it will be inter- 
esting to compare our results with more refined data in 
the future. 

This comparison can be interesting also for another 
reason. Indeed, the EPNJL model and the PNJLg can de- 
scribe the same QCD thermodynamics at zero and imag- 
inary quark chemical potential |3l| , but they differ qual- 
itatively for the interaction content. As we have shown, 
they have some quantitative discrepancy for what con- 
cerns the response to a strong magnetic field. Therefore, 
more refined Lattice data in magnetic field might help 
to discern which of the two models is a more faithful 
description of QCD. 

From our point of view, it is fair to admit that our 
study might have a weak point, namely, we miss a micro- 
scopic computation of the parameters ai, a2 in Eq. (O. 
For concreteness, we have used the best-fit values quoted 
in [sij . showing in addition that changing the value of 
To in the Polyakov loop effective potential as in [3l| , we 
obtain = Tl = 175 MeV, in excellent agreement with 
that reference. This is comforting, since it shows that 
our different UV-regulator does not affect the qualitative 
result at zero field, namely the coincidence of and T^. 
Besides, we are aware that Eq. ([5]) is just a particular 
choice of the functional dependence of the NJL coupling 
constant on the Polyakov loop expectation value. Differ- 
ent functional forms respecting the C and the extended 
Z3 symmetry are certainly possible, and without a rigor- 
ous derivation of Eq. ([5]) using functional renormalization 
group techniques as suggested in [s^l, it merits to study 
our problem using different choices for G{L) in the next 
future. For these reasons, we prefer to adopt a conserva- 
tive point of view: it is interesting that a model which 
is adjusted in order to reproduce Lattice data at zero 
and imaginary chemical potential, predicts that the two 
QCD transitions are entangled in a strong magnetic back- 
ground; however, this conclusion might not be definitive, 
since there exist other model calculations which share 
a common basis with ours, and which show a more pro- 
nounced split of the QCD transitions in a strong magnetic 
background. More refined Lattice data will certainly help 
to discern which of the two scenarios is the most favor- 
able. 

As a natural continuation of this work, it is worth to 
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perform the computation of the chiral magnetization |4l| 
at finite temperature. Besides, the technical machinery 
used here can be easily apphed to a microscopic study of 
the spectral properties of mesons in strong magnetic field. 
In this direction, it is interesting to compute the masses 
of the charged p-mesons at low temperature, in order to 
investigate their condensation at large magnetic field as 
suggested in psj . Furthermore, it would be interesting 
to make a complete (analytical or semi-analytical) study 
of the chiral limit, to estimate the effect of the magnetic 
field on the universality class of two-fiavor QCD. We will 
report on these topics in the next future. 
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